Unusual magnetic field-induced phase transition in the mixed state of 

superconducting NbSe2 



cn 

O 

o 



in 



o 

Oh' 



X3 

o 



> 

On 

m 

o 
m 
o 



I 

o 
o 



X 
S3 



A. V. Sologubenko^ I. L. Landau^-^, H. R. Ott\ A. Bilusic^-^, A. Smontara^, and H. Berger'^ 
^ Laboratorium fur Festkorperphysik, ETH Honggerherg, CH-8093 Zurich, Switzerland 
^Kapitza Institute for Physical Problems, 117334 Moscow, Russia 
^Institute of Physics, POB304, HR-10001 Zagreb, Croatia and 
"^Institut de physique de la matiere complexe, EPFL, CH-1015 Lausanne, Switzerland 

(Dated: February 2, 2008) 

The thermal conductivity k in the basal plane of single-crystalline hexagonal NbSe2 has been 
measured as a function of magnetic field H, oriented both along and perpendicular to the c-axis, at 
several temperatures below Tc. With the magnetic field in the basal plane and oriented parallel to 
the heat flux we observed, in fields well below Hc2, an unexpected hysteretic behavior of k,{H) with 
all the generic features of a first order phase transition. The transition is not manifest in the k{H) 
curves, if H is still in the basal plane but oriented perpendicularly to the heat-flux direction. The 
origin of the transition is not yet understood. 

PACS numbers; 74.60.Ec, 74.70.-b, 74.25.Fy 



The layered hexagonal compound NbSe2 has been of 
interest to researchers for some time, mainly because 
of the relatively strong anisotropy of its superconduct- 
ing properties. Pure samples exhibit a superconducting 
critical temperature Tc « 7.2 K, and the upper critical 
fields for H along and perpendicular to the hexagonal 
c-axis are H^^ = 40 kOe and = 140 kOe,^ respec- 
tively. The anisotropy of the electron system is also re- 
flected in the charge density wave (CDW) transition at 
TcDW = 34 K (in pure NbSe2), interesting by itself and 
also in view of the competition between the CDW for- 
mation and superconductivity.^ Another important fea- 
ture of NbSe2 is the existence of a variety of vortex mat- 
ter related phenomena, such as a pronounced peak effect 
and various phase transitions and instabilities (see, e.g., 
Refs. 3,4,5 and references therein). 

In this paper, we present the observation of an unex- 
pected hysteretic behavior of the thermal conductivity 
in external magnetic flelds k with the characteristic fea- 
tures of a first order phase transition. Several anomalies 
associated with the transition are yet to be explained. 

The single crystal used in the experiments was grown 
by employing the iodine vapor transport method as de- 
scribed in Ref. 6. Two bar-shaped samples with dimen- 
sions of 4.4 X 0.88 X 0.42 mm^ and 2.8 x 0.51 x 0.50 mm^, 
were cut from the crystal. Very detailed measurements 
of k(T, H) were done on one of the samples and the 
other specimen was used to ascertain that the observed 
effects described below are indeed sample-independent. 
The thermal conductivity was measured by employing 
the standard uniaxial heat flux method in an experimen- 
tal setup that was also used in our recent study of k(T, H) 
in MgB2.^'^ The heat flux Q was kept in the basal plane 
and along the longest lateral dimensions of the sample. 
The temperature difference AT along the specimen was 
kept to about 3 - 5 % of the average sample tempera- 
ture. The zero-field superconducting transition temper- 
ature Tc = 7.2 K was determined via four-contact dc re- 
sistivity measurements and from the results of measure- 
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FIG. 1: The basal plane thermal conductivity of single crys- 
talline NbSe2 as a function of differently oriented magnetic 
fields at T = 0.85 K. The open circles and squares correspond 
to increasing and decreasing fields, respectively. 



ments of the magnetic susceptibility, employing a com- 
mercial Quantum Design SQUID magnetometer. 

Fig. 1 shows k{H) aiT — 0.85 K for all three orienta- 
tions of the magnetic field that we have explored. Similar 
results were obtained at other temperatures in the inter- 
val between 0.38 K and Tc. The data for H \\ c are 
in good agreement with the recent results presented by 
Boaknin et al..^ The main features of k(T, H) can readily 
be understood by considering that the total thermal con- 
ductivity consists of two components, namely the phonon 
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contribution Kph and the quasiparticlc (electron) contri- 
bution Kg. The normal electronic excitations in the cores 
of the vortices in the mixed state considerably reduce the 
phonon mean free path, but they also add positively to 
the heat transport. Hence when the external magnetic 
field is enhanced to above the lower critical field Hci, 
Kph{H) is efficiently reduced but, simultaneously, Ke{H) 
starts to increase. The result is a minimum in the k{H) 
curve in a magnetic field slightly exceeding Hd- These 
features may be seen in Figs. 1 and 2(b). In the nor- 
mal state, above the upper critical field i?c2) the thermal 
conductivity varies only weakly with H . It reflects the 
field induced positive variation of the electrical resistiv- 
ity of NbSe2 in the normal state, which also affects the 
electronic contribution Ke to k. 

This general behavior of k{H) is observed for all ori- 
entations of the magnetic field, but, with H J- c and par- 
allel to the heat-fiux direction, a pronounced hysteretic 
behavior in k{H), as shown in Fig 2, is observed. This 
feature, emphasized in Fig. 2 (a), is a clear manifesta- 
tion of some kind of phase transition, most likely of first 
order. As we discuss in detail below, some details are 
not compatible with a conventional first order transition, 
however. Inspecting the data presented in Fig 2, it is 
clear that two different k{H) regimes have to be distin- 
guished. They correspond to two different modifications 
or phases of the mixed state, one of which is stable in 
low fields (^-phase) and the other (/i-phase) corresponds 
to the equilibrium state in high magnetic fields. In the 
following we shall use k'^''\H) and k^'^\H) to denote the 
thermal conductivities of the I- and /i-phases of the mixed 
state, respectively (sec Fig. 2(a) for the definitions of , 
K^'*' and the corresponding characteristic fields). While 
the ?-phasc is stable up to H = H2 , the high-field modifi- 
cation remains stable down to H = Hi < H^. Also note 
that the transition from one modification to the other is 
rather narrow upon field reduction, whereas the transi- 
tion in increasing field is more than twice as wide, with 
k{H) being almost constant in the transition region. All 
the features shown in Fig. 2 were confirmed to be re- 
producible and they were observed for both investigated 
samples. 

The hysteresis is completely absent if the magnetic 
field, while still parallel to the basal plane, is rotated to 
be perpendicular to the heat flux direction (see Fig 2(b)). 
This disappearance of the hysteresis is quite unexpected 
because, considering the hexagonal symmetry of the crys- 
tal structure of NbSe2 , no anisotropy is expected in the 
basal plane. It may also be seen that k{H) is notably 
higher if the magnetic field is oriented perpendicularly to 
the heat flux. 

Another unusual feature of this transition is the consid- 
erable width of the hysteresis loop with H2/H* w 1.7. It 

turns out that both H* and H2 increase with increasing 
temperature (see Fig. 3 (a)). Considering that all critical 
fields of a superconductor usually decrease with increas- 
ing temperature, this is again an unexpected observation. 
The magnetic field ranges H2 < H < H2 in increasing 
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FIG. 2: (a) Characteristic features of the hysteresis exempli- 
fied by k{H) at T = 0.38 K. The open circles and squares cor- 
respond to increasing and decreasing field, respectively, (b) 
The thermal conductivity as a function of magnetic field ori- 
ented perpendicularly to the c-axis at several temperatures. 
Note the difference for different field orientations. 



and < H < Hi in decreasing fields, respectively, 
correspond to transition regions where the I- and the h- 
phase coexist. Because both kM^{H) and k^'^''{H) may 
quite well be approximated by parallel straight lines, the 
volume ratio between the two phases for each point in 
the transition regions can easily be evaluated from the 
experimental data. 

Partial hysteresis loops, as shown in Fig. 4, are also un- 
usual. If a reduction of the magnetic field is started from 
the transition region between H2 and H2 (Fig. 4 (a)), 
the measured k,{H) curve is practically parallel to k^'^ (H) 
and K*^''^ (H) until Hi is reached. Roughly the same fea- 
ture emerges if the magnetic field is enhanced from the 
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FIG. 3: (a) Temperature dependencies for the fields and 
(b) Temperature dependence of the normahzed difference 
Ak = K^'^ - k'-''\ For definitions, see Fig. 2 (a). 
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FIG. 4: Partial hysteresis loops of k(H) a.t T = 0.85 K. The 
oijcii circles and squares represent the full loop and arc the 
same as in Fig. 2(b). (a) H is first increasing from to 12 kOe 
and then decreasing, (b) H is first enhanced from 30 kOe to 
8.8 kOe and then reduced. The broken vertical lines mark the 
special fields values. 



transition region Hi > H > H* (sec Fig. 4(b)). These 
observations imply that in the limited field range between 
the vertical broken lines in Fig. 4. the volume ratio be- 
tween the I- and the h-phase in the sample is independent 
of the applied magnetic field and can only be changed if 
the magnetic field exceeds H2 (in increasing field) or is 
reduced to below Hi (in decreasing field). This behavior 
is rather difficult to understand in terms of a common 
first order phase transition. Indeed, in conventional first 
order transitions, the main reason for a hysteresis is the 
inertia toward generating a sufficiently large nucleus of a 
new phase. Once the nucleus exists, the transition usu- 
ally happens in a rather narrow range of fields or tem- 
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FIG. 5: K vs time for two values of the external magnetic 
field at T = 0.85K. 



peratures. Here we have a completely different situation. 
Two different phases coexist in fixed amounts over an 
extended field region and it needs to be understood why 
the transformation from one phase to the other does not 
happen upon varying the applied magnetic field within 
the limits mentioned above. 

We also tested the stability of k at fixed values of the 
magnetic field in the transition regions and at constant 
temperature. As may be seen in Fig. 5, there are no 
observable changes of k with time. This observation pro- 
vides again clear evidence that the transformation of one 
phase into the other cannot happen as easily as it usually 
does in the case of common first order phase transitions. 
In an extended range of the magnetic field both the I- and 
the /i-phase are present and stable in the sample. Taking 
into account that near the edges of the hysteresis loop 
only one of them represents the equilibrium, the stabil- 
ity of the described situation implies that some potential 
barrier seems to prevent the completion of the transition 
from one equilibrium phase to the other. 

Before we discuss the situation further, we note that 
at temperatures as low as 0.38 K, the phonon contribu- 
tion Kph to the thermal conductivity in the mixed state is 
small compared to Kg, certainly in fields beyond the k{H) 
minima (see Figs. 1 and 2). Because wc are mainly inter- 
ested in low magnetic fields H < 0.1Hc2, no overlap of 
vortex cores needs to be considered and each vortex line 
is expected to contribute to the heat transport individ- 
ually. If the magnetic field is perpendicular to the heat 
fiux, the vortex motion along the heat fiux, induced by a 
temperature gradient VT, also contributes to k, but this 
contribution is not expected to be significant (see Ref. 10 
and references therein). 

The results presented above clearly imply the existence 
of a phase transition in the mixed state of our NbSe2 
samples with the following main features: 

(i) A very large difference Ak = k^'^ — k^'^^ between 
the thermal conductivities in the low and the high field 
phase. If we assume that Ak is only due to a change of 
the vortex density, we have to assume an unrealistically 
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large variation of the magnetic induction AB w 3 kG 
(see Fig. 2 (a)). 

(ii) A very large width of the hysteresis loop with 
H*/H* = 1.7. 

(iii) The coexistence of the I- and the /i-phase of the 
mixed state in an extended range of magnetic fields. The 
adopted volume ratio is practically independent of H in 
this range (Fig. 3). This particular feature is incompati- 
ble with a conventional first order transition. 

(iv) A completely different behavior of k{H), if the 
magnetic field is oriented perpendicularly to the heat 
flux, exhibiting no visible signs of a phase transition 
(Fig. 2 (b)). 

(v) A shift of the transition to higher fields with in- 
creasing temperature (Figs. 2 (b) and 3 (a)). 

Although the phase transition manifests itself by a 
rather prominent feature in the k{H) curves, we could 
not find a reasonable scenario which would explain all 
the different features of the transition. Below we briefiy 
consider several possible phase transitions in the mixed 
state of type-II superconductors and discuss their rele- 
vance for our experimental observations. 

1. The well known vortex- lattice melting transition ex- 
plains neither the magnitude of the observed difference 
Ak = k'-'^ — K^'^\ nor the vanishing of the transition for 
if _L Q. Quite unlikely are other vortex ordering tran- 
sitions considered in the literature, because all of them 
invoke only small changes of the vortex density. We also 
recall that the vortex lattice melting is a typical first or- 
der transition and therefore it is incompatible with the 
partial hysteresis loops shown in Fig. 4. 

2. It was pointed out that in superconductors with 



other ordering phenomena, a first order phase transi- 
tion may occur due to the competition between super- 
conductivity and the additional ordering phenomenon. 
As an example of additional ordering, the authors of 
Ref. 11 considered a charge density wave formation. This 
matches exactly the case for NbSe2 but, because the as- 
sumed transition is of first order, it is not clear how the 
coexistence of two mixed state phases and the depen- 
dence of the k{H) feature on the orientation of the mag- 
netic field can be explained. 

3. The first order transition discussed in Ref. 5 was 
observed for -ff || c and therefore seems not to be relevant 
to our case. The same applies to other reports of insta- 
bilities of the mixed state of NbSe2 which were observed 
either close to Hc2 or for other field orientations.^''' 

4. The possible phase transition discussed in Ref. 12 is 
expected to occur in fields rather close to Hc2 and there- 
fore is a priori considered as unlikely to be the reason for 
our observation. 

In conclusion, measurements of the thermal conductiv- 
ity in external magnetic fields revealed an unusual phase 
transition in the mixed state of NbSe2. The origin and 
the nature of this transition is not yet understood and 
further experiments arc needed to clarify the situation. 
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